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Abstract 

Internal jet structure in dijet production in deep-inelastic scattering is measured with the 
HI detector at HERA. Jets with transverse energies £r,Breit > 5GeV are selected in the 
Breit frame employing k± and cone jet algorithms. In the kinematic region of squared 
momentum transfers 10 < Q 2 < 120 GeV 2 and x-Bjorken values 2 • 10 -4 < XBj < 
8 • 10~ 3 , jet shapes and subjet multiplicities are measured as a function of a resolution 
parameter. Distributions of both observables are corrected for detector effects and presented 
as functions of the transverse jet energy and jet pseudo-rapidity. Dependences of the jet 
shape and the average number of subjets on the transverse energy and the pseudo-rapidity 
of the jet are observed. With increasing transverse jet energies and decreasing pseudo- 
rapidities, i.e. towards the photon hemisphere, the jets are more collimated. QCD models 
give a fair description of the data. 
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1 Introduction 



A sizeable fraction of the final states produced in high energy collisions shows the characteristic 
feature of large amounts of hadronic energy in small angular regions. These collimated sprays 
of hadrons (called jets) are the observable signals of underlying short distance processes and 
are considered to be the footprints of the underlying partonic final states. Quantitative studies 
of jet production require a precise jet definition, which is given by a jet finding algorithm. 

Jets so defined exhibit an internal structure which is sensitive to the mechanism by which a 
complex aggregate of observable hadrons evolves from a hard process. The understanding of 
this mechanism involves higher orders of the strong coupling constant in perturbation theory as 
well as non-perturbative contributions. This is a challenging task for theory. Recently, for some 
specific hadronic final state quantities, encouraging results have been obtained by exploiting 
the characteristic power behaviour of non-perturbative effects and by analytical, approximate 
calculations of perturbative QCD parton evolution down to the semi- soft regime Fur- 
thermore, since jet production rates are used to test the predictions of perturbative QCD, the 
understanding of their detailed properties and internal structure is an important prerequisite. 

The internal structure of jets has been studied in e + e [|3|] and in hadron-hadron colli- 
sions [Q], At the e ± p collider HERA, these investigations can be performed in photoproduction 
(Q 2 ?s GeV 2 ) and in deep-inelastic scattering (DIS) at large squared four momentum trans- 
fers Q 2 . In a previous publication we have measured the E T dependence of the jet width [|5|] 
in photoproduction. Recently, the ZEUS collaboration has investigated jet shapes in photopro- 
duction [§] and in DIS at Q 2 > 100 GeV 2 [JTp. Both analyses are carried out in the laboratory 
frame. This means that for DIS at high Q 2 mostly events with only one jet enter the analysis. 

The hadronization of the current jet in deep-inelastic scattering in the Breit frame has al- 
ready been studied with event shape variables [^p, charged particle multiplicities and fragmen- 
tation functions [0]. In this paper we take the first steps towards a complete understanding of 
jet properties in DIS. We analyse the hadronization of jets in multijet production in the Breit 
frame. The Breit frame, where the virtual photon interacts head-on with the proton, has been 
chosen in this analysis because here the produced transverse)] energy, -£/T,Breitj directly reflects 
the hardness of the underlying QCD process. We present measurements of internal jet struc- 
ture in a sample of inclusive dijet events with transverse jet energies of -Er.Breit > 5 GeV, 
10 < Q 2 < 120 GeV 2 and 2 ■ 10~ 4 < x Bj < 8 ■ 10~ 3 . This is the £ T ,Breit range where jet cross 
section measurements are currently performed at HERA and compared to perturbative QCD 
calculations (e.g. [[TTH, [TT]]). The analysis is based on data taken in 1994 with the HI detector 
at HERA when 27.5 GeV positrons collided with 820 GeV protons. The data correspond to an 
integrated luminosity of C- in t — 2pb _1 . 

Jets are defined in the Breit frame by kj_ and cone jet algorithms. Two observables, jet 
shapes and, for the first time, subjet multiplicities, are studied. The jet shape measures the 
radial distribution of the transverse jet energy around the jet axis. For the k± cluster algorithm 
we have also measured the multiplicity of subjets, resolved at a resolution scale which is a 

'transverse with respect to the z-axis which is given by the axis of the virtual photon and the proton. 
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fraction of the jet's transverse energy. Both observables are presented for different ranges of the 
transverse jet energy and the pseudo-rapidity^ of the jets in the Breit frame. 

The paper is organized as follows. Section ^ gives a brief description of the HI detector. 
In section ^| we introduce the jet algorithms used in the analysis and give the definition of the 
measured observables in section |j. In section ^| we give a short description of the QCD models 
which are used for the correction of the data and to which the results are later compared (in 
section The data selection and the correction procedure are described in sections |^ and [7| 
and the results are discussed in section ||. 



2 The HI Detector 

A detailed description of the HI detector can be found elsewhere flT2|]. Here we briefly introduce 
the detector components relevant for this analysis: the liquid argon (LAr) calorimeter [|13|], the 
backward lead-scintillator calorimeter (BEMC) [fl^Q, and the tracking chamber system [15]. 

The hadronic energy flow is mainly measured by the LAr calorimeter extending over the 
polar angular range 4.4° < 9 < 154° with full azimuthal coverage. The polar angle 9 is 
defined with respect to the proton beam direction (+z axis). The LAr calorimeter consists of an 
electromagnetic section (20 — 30 radiation lengths) with lead absorbers and a hadronic section 
with steel absorbers. The total depth of both calorimeters varies between 4.5 and 8 interaction 
lengths. Test beam measurements of the LAr calorimeter modules show an energy resolution 



of a E /E « 0.50/a/-E [GeV] © 0.02 for charged pions \\L($. The absolute scale of the hadronic 
energy is known for the present data sample to 4%. 

The scattered positron is detected by the BEMC with a depth of 22.5 radiation lengths 
covering the backward region of the detector, 155° < 9 < 176°. The electromagnetic energy 
scale is known to an accuracy of 1%. 

The calorimeters are surrounded by a superconducting solenoid providing a uniform mag- 
netic field of 1.15 T parallel to the beam axis in the tracking region. 

Charged particle tracks are measured in two concentric jet drift chamber modules (CJC), 
covering the polar angular range 15° < 9 < 165°. The forward tracking detector covers 7° < 
9 < 25° and consists of drift chambers with alternating planes of parallel wires and others 
with wires in the radial direction. A backward proportional chamber (BPC) with an angular 
acceptance of 151° < 9 < 174.5° improves the identification of the scattered positron. The 
spatial resolution for reconstructed BPC hits is about 1 .5 mm in the plane perpendicular to the 
beam axis. 



3 Jet Definitions 

The jet algorithms used in this analysis are applied to the particles boosted into the Breit frame. 
Particle refers here either to an energy deposit in the detector (see section ^|), to a stable hadron 

2 The pseudo-rapidity r\ is defined as r\ = — ln(tan 8/2) where 9 is the polar angle with respect to the proton 
direction. This definition is chosen in both the laboratory frame and the Breit frame. 
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or a partem in a QCD model calculation. In all cases the scattered positron is excluded. The 
Breit frame is defined by q + 2x-q- s P = 0, where q and P are the momenta of the exchanged 
boson and the incoming proton. The z-axis is defined as the direction of the incoming proton. 

In the following analysis we use two different jet definitions: a cone algorithm and a k± 
cluster algorithm. Both jet definitions are invariant under boosts along the z-direction. The 
recombination of particles is carried out in the E T recombination scheme, which is based on 
transverse energies E T , pseudo-rapidities 77 and azimuthal angles <\> of the particles. The trans- 
verse energy and the direction of a jet are defined by 

E T jet = ^ E T ,j, Tftet = p' - , 0jet = ^ ~ , (1) 

where the sums run over all particles i assigned to the jefl 



3.1 Cone Algorithm 

Based on the original proposal of Sterman and Weinberg [Jl7|] many different implementations of 
cone algorithms have been developed. While the basic idea of the cone algorithm is simple and 
very intuitive, an operational definition is non-trivial. The resulting jet cross sections depend on 
how the algorithm treats the choice of jet initiators and configurations of overlapping jet cones. 
It has repeatedly been pointed out that many definitions of cone algorithms are not infrared 
and/or collinear safe [[T^,[T9[]. 



In this analysis we use the definition implemented in the algorithm PXCONE [ |20| ] which 
does not suffer from the problems discussed in [P"8||I9[]. This definition, which corresponds 
closely to the Snowmass proposal [|21|] and to the algorithm used in the CDF experiment 



is also used by the OPAL collaboration Q23Q. 



Particles are assigned to jets based on their spatial distance R in pseudo-rapidity and azimuth 
space (R 2 = Ar/ 2 + A</> 2 ). The algorithm operates as follows: 



1. Each particle is considered as a seed of a jet, for which steps 2-4 are performed. 

2. The jet momentum is calculated from all particles within a cone of radius R around the 
seed direction using eq. (jl]). 

3. If the jet direction differs from the seed direction, the jet direction is taken as the new 
seed direction and step 2 is repeated. 

4. When the jet direction is stable the jet is stored in the list of "protojets" (if it is not identical 
with a protojet already found). 

5. The steps 2 to 4 are repeated for all midpoints of pairs of protojets as seed directions^. 
This leads to the infrared safety of the procedure JT^]. 

3 All particles are considered massless by setting Ei — \pl\. 

4 In practice it is sufficient to do this only for pairs of protojets with a distance between Rq and 2 Rq. 
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6. Protojets with transverse energies of E T ^ ct < e are removed from the list. The cut-off 
parameter e specifies below which transverse energies protojets are not considered in the 
overlap treatment (steps 7-8). 

7. All remaining protojets that have more than a fraction / of their transverse energy con- 
tained in a protojet of higher transverse energy are deleted. 

8. All particles that are contained in more than one protojet are assigned to the protojet 
whose center is nearest in (77, <p). 

9. The jet momenta are recalculated using eq. ([[]). All protojets with E T ^ et < e are deleted 
and the remaining ones are called jets. 

The jets with the highest transverse energies are considered in the analysis. Due to the 
reassignment of particles to jets and the recalculation of the jet axis (steps 7, 8) it may happen 
that single particles within a jet have a distance larger than _R t0 the jet axis. This analysis is 
made with the parameter settings e = 5 GeV, / = 0.75 and a cone radius of R = 1.0. 



3.2 Inclusive k± Algorithm 

The ambiguities that occur for cone jet definitions (choice of seeds, overlapping cones) are 
avoided in cluster algorithms which successively recombine particles to jets. One definition 
of such an algorithm (proposed in p4) ] and implemented in the KTCLUS algorithm []25|]) has 
properties very similar to cone algorithms. As in the cone algorithm the clustering procedure 
is based on the longitudinally boost-invariant quantities E T , Ar], A0. The minimum of all dis- 
tances between particles is determined and either the corresponding pairs of particles are merged 
into pseudo-particles or single (pseudo-) particles are declared as jets. This process is iterated 
until no particles are left: 

1. We start with a list of all particles and an empty list of jets. 

2. For each particle i as well as for each pair of particles the distances di and dij are 
calculated 

d i = El l R 2 and c% = min^, E%J) R% with R% = A77J + &<t>% . (2) 

3. The smallest value of all the di and d^ is labeled d min . 

4. If rf min belongs to the set of d^, the particles i and j are merged into a new particle using 
the recombination prescription in eq. ([]]) and removed from the list of particles. 

5. If d mhl belongs to the set of di, the particle i is removed from the list of particles and 
added to the list of jets. 

6. When no particles are left (i.e. all particles are included in jets) the procedure is finished. 

The last jets that entered the list are the ones with highest transverse energies. These jets 
are considered in the analysis. This jet definition implies that particles with Rij < R are 
subsequently merged, so that all final jets are separated by distances R^ > Rq. It is still 
possible that particles inside a jet have a distance R^ > Ro to the jet axis and that particles with 
Rij < R are not part of the jet. The parameter _R is set to R = 1.0. 
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4 The Observables 



Two observables of internal jet structure are investigated in this analysis. They are sensitive to 
different aspects of jet broadening. 

The jet shapes are studied for the cone and the k± algorithm. This observable measures the 
radial distribution of the transverse jet energy only and is affected by hard and by soft processes 
over the whole radial range. 

A natural choice for studying the internal structure of jets with the k± cluster algorithm is 
the multiplicity of subjets, resolved at a resolution scale which is a fraction of the jet's transverse 
energy. These subjet multiplicities are sensitive to more local structures of relative transverse 
momentum within a jet. Here the perturbative and the non-perturbative contributions are better 
separated. While at larger values of the resolution parameter perturbative contributions domi- 
nate, at smaller values non-perturbative contributions become increasingly important. 



The jet shape \&(r) is defined as the fractional transverse jet energy contained in a subcone of 
radius r concentric with the jet axis, averaged over all considered jets in the event sample 



where iVj Cts is the total number of these jets. As proposed in [|TP|], only particles assigned by the 
jet algorithm to the jet are considered. 

Usually the denominator in the definition of ^ is given by the summed E T of all particles 
within a radius R to the jet axis. This means that ^(r/R = 1) = 1. In our definition @ of 
^ the denominator is given by the transverse energy of the jet. Since neither for the cone nor 
for the k± definition are all particles necessarily assigned to a jet within a radius of r/R < 1 
to the jet axis, ^(r/R = 1) is not constrained to have the value of one. With this choice of our 
observable we are also sensitive to the amount of transverse jet energy outside the radius Ro. 

4.2 Subjet Multiplicities 

For each jet in the sample the clustering procedure is repeated for all particles assigned to the 
jet. The clustering is stopped when the distances yij between all particles i, j are above some 
cut-off y cut 



and the remaining (pseudo-)particles are called subjets. The parameter y cut defines the minimal 
relative transverse energy between subjets inside the jet and thus determines the extent to which 
the internal jet structure is resolved. From this definition it follows that for y cut > 0.25 no subjet 
is resolved (therefore the number of subjets is one), while for y cut — > every particle in the jet 
is a subjet. The observable that is studied in this analysis is the average number of subjets for a 
given value of the resolution parameter, for values y cut > 1CT 3 . 



4.1 The Jet Shape 




(3) 



mmjE^E^) (A^ + Aj|) 



> Vent 



(4) 
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5 QCD Models 



A simulation of the detailed properties of the hadronic final state is available in the form of 
Monte Carlo event generators. They include the matrix element of the hard subprocess in first 
order of the strong coupling constant a s , approximations of higher order QCD radiation effects, 
and a model to describe the non-perturbative transition from partons to hadrons. 

The LEPTO Monte Carlo [^5p incorporates the 0(a s ) QCD matrix element and takes higher 
order parton emissions to all orders in a s approximately into account using the concept of parton 
showers [ jZ7| ] based on the leading logarithm DGLAP equations [^Sp. QCD radiation can occur 
before and after the hard subprocess. The formation of hadrons is performed using the LUND 
string model implemented in JETSET [TJ(J. 

The HERWIG Monte Carlo [ [31"1 ] also includes the 0(a s ) QCD matrix element, but uses 
another implementation of the parton shower cascade which takes coherence effects fully into 
account. The hadronization is simulated with the cluster fragmentation model 



In ARIADNE [ TJ3| ] gluon emissions are treated by the colour dipole model [TMp assuming a 
chain of independently radiating dipoles spanned by colour connected partons. The first emis- 



sion in the cascade is corrected to reproduce the matrix element to first order in a s Q35Q. 



DJANGO [p6fl provides an interface between the event generators LEPTO or ARIADNE 



and HERACLES Q37Q which makes it possible to include 0{a) QED corrections at the lepton 
line. 



6 Data Selection 



The analysis is based on HI data taken in 1994 corresponding to an integrated luminosity of 
d n t — 2pb _1 . The event selection closely follows that described in a previous publication [[T0|]. 
DIS events are selected where the scattered positron is measured in the acceptance region of 
the BEMC at energies where trigger efficiencies are approximately 100%. To ensure a good 
identification of the scattered positron and to suppress background from misidentified photo- 
production events the following cuts are applied: 



• The cluster of the positron candidate must have an energy-weighted mean transverse ra- 
dius below 5 cm. 

• A reconstructed BPC hit within 5 cm of the straight line connecting the shower center 
with the event vertex is required. 

• The z position of the reconstructed event vertex must be within ±30 cm of the nominal 
position. 

• A cut on 35 GeV < ^2(E — p z ) < 70 GeV is applied, where the sum runs over all energy 
deposits in the calorimeter. In neutral current DIS events without undetected photon 
radiation the quantity ^2(E — p z ) is expected to be equal to twice the energy of the initial 
state positron. This cut reduces the contribution from photoproduction events as well as 
events where hard photons are radiated collinear to the incoming positron. 
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The event kinematics are calculated from the polar angle 9 e i and the energy E' el of the 
scattered positron via Q 2 el = 2E E' el (l + cos0 e j)> Vei — 1 — E' el /(2Eo){l — cos9 e {) and 
^Bj = Q 2 /(sy). Eq denotes the energy of the incoming positron and s the ep centre-of- 
mass energy squared. Events are only accepted, if E' el > 11 GeV, 156° < 9 e i < 173°, 
Q 2 > 10 GeV 2 and y > 0.15. The resulting kinematic range is 10 < Q 2 < 120 GeV 2 and 
2 ■ 10~ 4 < x Bi < 8 ■ 10- 3 . 

Jets are defined by the algorithms described in section |3|. The input for the jet algorithms 
consists of a combination of energy clusters from the calorimeter and track momenta measured 
in the central and forward trackers (as described in [[TO]]). While all energy clusters are consid- 
ered, the four momentum of each single track is only allowed to contribute up to a momentum 
of 350 MeV. This procedure partly compensates for energy losses in the calorimeter due to 
dead material and noise thresholds. It reduces the dependence of the jet finding efficiency on 
the pseudo-rapidity of the jet and improves the reconstruction of the transverse jet energy p8|]. 

The objects from tracking and calorimeter information are boosted to the Breit frame where 
the jet algorithms are applied. We select events with at least two identified jets with transverse 
energies of £r,Breit > 5 GeV in —1 < 77jet,iab < 2. The two jets with the highest E T ,Bvdt are 
considered in the analysis. The event sample for the inclusive k± algorithm (the cone algorithm) 
consists of 2045 (2657) dijet events. 



7 Correction of the Data 



The data are corrected for detector effects and QED radiation from the lepton. The detector 
response is determined using events from Monte Carlo event generators that were subjected to 
a detailed simulation of the HI detector. The following event generators are used: ARIADNE 
interfaced in DJANGO (with and without the inclusion of QED corrections) and LEPTO. Both 
generators give a good description of the kinematic variables of the inclusive DIS data sample 



as well as of the angular and transverse energy distributions of the jets 0390. We also observe a 



reasonable description of the observables introduced in section |] (see section ^[). 

The measured data points are corrected bin-by-bin for detector effects. Using the generated 
event samples, the correction factor for each bin is determined as the ratio of the generated value 
of the observable and the value that is reconstructed after detector simulation. These correction 
factors are independent of the inclusion of QED radiation effects as included in DJANGO. 
Their dependence on details of the modeling of the hadronic final state is taken into account 
by considering the difference between the correction factors from ARIADNE and LEPTO as 
systematic uncertainty. 

For the k± (cone) algorithm the corrections for \I/(r) are below 10 % (13 %) for subcone radii 
r > 0.3 and always below 27% (23%). The corrections for (N snb j et (y cut )) are always below 
7%. The correction factors from both QCD models are in good agreement (they differ typically 
by not more than 2 %) for the jet shapes as well as for the subjet multiplicities The final 
correction factors are taken to be the mean values of the two models, taking the spread as the 
error. In addition we have varied the calibration of the hadronic energy scale in the data sample 
in the range of ±4% around the nominal value. The error is estimated as the maximal deviation 
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from the results at the nominal value. For all observables it is at most 5%. The overall systematic 
error is calculated by adding the errors from the model dependence and from the uncertainty of 
the hadronic energy scale in quadrature. In all figures the statistical and systematic errors are 
added in quadrature. Since each jet enters in all bins of a distribution, all errors are correlated. 

The background from misidentified photoproduction events is estimated with a sample of 
photoproduction events generated with PHOJET [ |TO| ] and is found to be negligible. 

8 Results 

The jet shape and the subjet multiplicity are presented as functions of quantities directly related 
to the single jets, namely the transverse jet energy (£/r,B re it) and the pseudo-rapidity (r] B reit) in 
the Breit frame. We also investigated whether the observables depend on the event kinematics. 
The jet shapes and subjet multiplicities were compared for two bins of Q 2 (Q 2 < 20 GeV 2 and 
Q 2 > 20 GeV 2 ) and xbj (£ B j < 8 • 10~ 4 and x B j > 8 • 10~ 4 ) respectively. No dependence on 
Q 2 and x B j has been observed. 

8.1 Jet Shapes 

The radial dependence of the jet shape \&(r) for the kj_ algorithm is shown in Fig. [I] in different 
ranges of the pseudo-rapidity in the Breit frame. The results for jets of transverse energies 
5 < -&r,Breit < 8 GeV and -Sr.Breit > 8 GeV are superimposed. The jet shape ^(r) increases 
faster with r for jets at larger transverse energies, indicating that these jets are more collimated. 
The same tendency is seen for the jets defined by the cone algorithm which are compared to the 
jets found by the k± algorithm in Fig. Q For both jet definitions we also observe a dependence 
of the jet shape on the pseudo-rapidity of the jets. Jets towards the proton direction (at larger 
values of r/Breit) are broader than jets towards the photon direction (smaller rfereit)- In the region 
where the jets are most collimated (-&r,Breit > 8 GeV and r/Breit < 2.2), very similar jet shapes 
are observed for the k± and cone algorithms. The broadening of the jets for smaller i?T,Breit and 
larger r/Breit is more pronounced for the cone jet definition. 

Recently jet shapes have been measured in dijet production in photon-photon collisions [ pTT| ] 
for jets defined by a cone algorithm at transverse energies comparable to those presented here. 
The jet shapes in photon-photon collisions (where no t] dependence is observed) are very similar 
to those measured in DIS in the Breit frame at 77 Bre it < 1.5. 

8.2 Subjet Multiplicities 

The subjet multiplicities for the k± algorithm are displayed in Fig. [3|. The average number of 
subjets (iV su bjet(|/cut)) as a function of the subjet resolution parameter at y cut > 10~ 3 is plotted. 
Towards smaller values of y cut , an increasing number of jet fragments with smaller relative 
transverse momenta is resolved. The number of subjets at a given value of y cut reflects the 
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amount of relative transverse momentum with respect to the jet axis. The subjet multiplicity is 
therefore a measure of the broadness of the jet. 

At y cut = 1CT 3 a jet is on average resolved into 4.1 - 4.6 subjets, depending on -Er.Breit an d 
^Breit of the jelQ. For almost all values of y cut the subjet multiplicity is larger for jets at smaller 
-Er.Breit and larger r] Bieit , indicating broader jets. 

A summary of the results for both observables is given in Fig. 0. Here the -Er.Breit and the 
77Breit dependence of the jet shape and the average number of subjets are shown at an intermediate 
value of the resolution parameter (jet shape: r = 0.5 and subjet multiplicity: y cut = 10~ 2 ). 

Although the subjet multiplicities are sensitive to the jet broadening in a different way than 
the jet shapes, consistent conclusions can be drawn for both measurements. The jet broadening 
depends on both the transverse jet energy as well as the pseudo-rapidity in the Breit frame. 
While the pseudo-rapidity dependence is most pronounced at smaller transverse jet energy, the 
transverse energy dependence is stronger in the forward region (at larger pseudo-rapidities). 



9 Comparison with QCD Model Predictions 

The predictions of different QCD models are compared in Fig. |5] to the jet shapes measured for 
the k ± algorithm. The models LEPTO, ARIADNE and HERWIG all show E T)Bieit and r/ Breit de- 
pendences similar to that seen in the data. LEPTO gives the best description of the measured 
shapes for ?7 Breit < 2.2 while at ?7 Breit > 2.2 the predicted jet shapes are too broad. A reasonable 
description is also obtained by the ARIADNE model except for jets at smaller pseudo-rapidities 
where the jet shapes have the tendency to be too narrow. For the HERWIG model the jet shapes 
are narrower than those in the data in all E TjBreit and ^ B reit regions. The same observations as 
above are made when comparing these QCD models with the subjet multiplicities and with the 
jet shapes for the cone algorithm (not shown here). 

In QCD models the evolution of a jet is described by perturbative contributions (radiation of 
partons) and non-perturbative contributions (hadronization). Studies based on the LEPTO and 
HERWIG parton shower models show that all observables studied in this analysis are strongly 
influenced by hadronization. This process has the largest impact on the jet broadening in our 
kinematic region (Fig. ^). Basic characteristics of the perturbative contributions are however 
still visible after hadronization. The model prediction suggests that the large difference between 
quark and gluon-initiated jets before hadronization survives the hadronization process. This 



especially applies to jets with large transverse energies [391. 



Fig. ^| shows the jet shapes and the subjet multiplicities as predicted by the LEPTO parton 
shower model for the kj_ algorithm, separately for quark and gluon jets at £r,Breit 

> 8GeV 

and ?7Breit < L5. Gluon jets are broader than quark jets. The same prediction is obtained by 
the HERWIG parton shower model. Although the jets in HERWIG are slightly narrower, the 
differences between gluon and quark jets are equally large. In the phase space considered here, 
LEPTO and HERWIG (in agreement with next-to-leading order calculations) predict a fraction 

5 On average the jets in the data (as in the simulated events) consist of eleven calorimetric energy clusters. For 
the LEPTO generator this is also approximately the average multiplicity of stable particles inside the jets. 
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of approximately 80% photon-gluon fusion events with two quarks in the partonic final state. 
The jet samples of these models are therefore dominated by quark jets. Both model predictions 
for the jet shapes and the subjet multiplicities therefore mainly reflect the properties of the quark 
jets as can be seen in Fig. ^. These predictions give a reasonable description of the data. Thus, 
we conclude, that the jets we observe are consistent with being mainly initiated by quarks. 



10 Summary 

Measurements of internal jet structure in dijet events in deep-inelastic scattering in the kinematic 
domain 10 < Q 2 < 120 GeV 2 and 2 ■ 10~ 4 < s B j < 8 • 10 -3 have been presented. Jet shapes and 
subjet multiplicities have been studied for jets of transverse energies -&r,Breit > 5 GeV defined 
by k±_ and cone jet algorithms in the Breit frame. 

The radial dependence of the jet shape and the dependence of the average number of subjets 
on the subjet resolution parameter y cut are both sensitive to different aspects of jet broadening. 
For both observables a dependence of the jet broadness on the transverse energy -Er,Breit and on 
the pseudo-rapidity in the Breit frame r] Bieit is seen. With increasing E TtBmit jets are narrower. 
Jets of the same -Er,Breit become broader towards the proton direction. This effect is more 
pronounced at lower .Sr.Breit- 

At lower -&r,Breit jets defined by the k± algorithm are more collimated than jets defined by 
the cone algorithm, while at higher -Er,Breit both algorithms produce very similar jets. 

The QCD models LEPTO, ARIADNE and HERWIG roughly reproduce the dependence 
of the jet shape and the subjet multiplicities on -Er.Breit and ^ Breit as seen in the data. LEPTO 
has a tendency to produce broader jets in the proton direction than measured. HERWIG and 
ARIADNE produce jets which are too collimated especially at higher transverse energies. We 
have reported earlier that in the same kinematic domain the predicted jet rates from LEPTO and 
HERWIG are about a factor of two below the data [ |10| ] . Since these models are able to reproduce 
the internal jet structure, this failure must be largely connected to an inadequate modeling of 
the underlying hard partonic subprocess. 

According to the parton shower models LEPTO and HERWIG, quark and gluon initiated 
jets differ both at the parton and at the hadron level. Both models predict that the jet sample is 
dominated by quark initiated jets. Since these models describe our data, we conclude that the 
observed jet structures are compatible with those of quark initiated jets. 
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Figure 1: The dependence of the jet shapes on the transverse jet energy in three pseudo-rapidity 
regions. The jet shapes at higher and at lower transverse jet energies for the inclusive k± algo- 
rithm are overlaid. The comparison is shown as a function of the jet pseudo-rapidity in the Breit 
frame (positive pseudo-rapidities are towards the proton direction). 
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Figure 2: Comparison of the jet shapes for the inclusive k± algorithm and the cone algorithm. 
The data are shown as a function of the transverse jet energy and the jet pseudo-rapidity in the 
Breit frame (positive pseudo-rapidities are towards the proton direction). 
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Figure 3: The average number of subjets as a function of the resolution parameter y cut for the 
inclusive k± algorithm. The data are shown in r^reit bins for different E T ,Bmit- 
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Figure 4: The jet shape \J/(r) at a fixed subcone radius r = 0.5 in three ry Br eit ranges and two 
-Er,Breit bins for the inclusive k± algorithm (left). The average subjet multiplicity at a fixed 
resolution parameter y cut = 0.01 in three r? B reit ranges and two -Er,Breit bins for the inclusive k± 
algorithm (right). Positive pseudo-rapidities are towards the proton direction. 
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Figure 5: The jet shapes for the inclusive k± algorithm. The data are shown as a function of the 
transverse jet energy and the jet pseudo-rapidity in the Breit frame (positive pseudo-rapidities 
are towards the proton direction). The results are compared to predictions of QCD models. 
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Figure 6: Model predictions of the internal structure of quark and gluon jets for the inclusive k± 
algorithm by the LEPTO parton shower model. The jet shapes (left) and the subjet multiplicities 
(right) are shown separately for quark and gluon induced jets with £r,Breit > 8 GeV and r?Brcit < 
1.5, together with the sum of both and the comparison to the HI measurement. The distributions 
of the observables before hadronization are also shown. 
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